home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / zdrot.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.3 KB  |  52 lines
  1.       SUBROUTINE ZDROT( N, ZX, INCX, ZY, INCY, C, S )
  2. *
  3. *     applies a plane rotation, where the cos and sin (c and s) are
  4. *     double precision and the vectors zx and zy are double complex.
  5. *     jack dongarra, linpack, 3/11/78.
  6. *
  7. *     .. Scalar Arguments ..
  8.       INTEGER           INCX, INCY, N
  9.       DOUBLE PRECISION  C, S
  10. *     ..
  11. *     .. Array Arguments ..
  12.       COMPLEX*16        ZX( 1 ), ZY( 1 )
  13. *     ..
  14. *     .. Local Scalars ..
  15.       INTEGER           I, IX, IY
  16.       COMPLEX*16        ZTEMP
  17. *     ..
  18. *     .. Executable Statements ..
  19. *
  20.       IF( N.LE.0 )
  21.      $   RETURN
  22.       IF( INCX.EQ.1 .AND. INCY.EQ.1 )
  23.      $   GO TO 20
  24. *
  25. *        code for unequal increments or equal increments
  26. *          not equal to 1
  27. *
  28.       IX = 1
  29.       IY = 1
  30.       IF( INCX.LT.0 )
  31.      $   IX = ( -N+1 )*INCX + 1
  32.       IF( INCY.LT.0 )
  33.      $   IY = ( -N+1 )*INCY + 1
  34.       DO 10 I = 1, N
  35.          ZTEMP = C*ZX( IX ) + S*ZY( IY )
  36.          ZY( IY ) = C*ZY( IY ) - S*ZX( IX )
  37.          ZX( IX ) = ZTEMP
  38.          IX = IX + INCX
  39.          IY = IY + INCY
  40.    10 CONTINUE
  41.       RETURN
  42. *
  43. *        code for both increments equal to 1
  44. *
  45.    20 DO 30 I = 1, N
  46.          ZTEMP = C*ZX( I ) + S*ZY( I )
  47.          ZY( I ) = C*ZY( I ) - S*ZX( I )
  48.          ZX( I ) = ZTEMP
  49.    30 CONTINUE
  50.       RETURN
  51.       END
  52.